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ABSTRACT 


Conventional nuclear reactors suffer from many drawbacks such 
as cofiplex reactor desi^ resulting in frequent safety problems, 
overall high cost of power generation sue to long construction 
times and retrofittings to comply with safety stipulations, etc. 
To overcome these problems, recently there has been considerable 
interest to develop an improved nuclear reactor. The Boiling Water 
Reactor using natural circulation as a mode for coolant flow is 
one of such new reactor vrtiich has the potential for overcoming 
many safety related problems due to elimination of external punps 
and simplify the reactor vessel design considerably. 

The present thesis analyses the feasibility of the natural 
circulation BWR for hi^r power rates satisfying one of the 
thermal hydraulic requirements, that is the core stability. The 
study shows that the natural circulation BWR for powers >1000 Mrf 


is indeed feasible. 
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SYMBOLS 

M : Mass of the fuel rod (Kg), 

p : Operating pressure of reactor (MPa). 

: Specific heat of fuel tnaterlad (KJ/Kg/ °C). 

: Fuel temperature (*C). 

: Coolant temperature (^^C). 

: Steady state fuel temperature (°C). 

: Steady state coolant temperature (*C). 

: Perturbation In fuel temperature (°C). 

: Perturbation in coolant temperature (*C) 

q : Fuel power d^slty (KW/1). 

q" s Heat flux (MW/m^). 

3 

V : Volume of fuel rod (m ). 

A : Surface area of fuel rod (m^). 

q^ : Steady state fuel power density (WV/l). 

q“ : Steady state heat flux (MW/in^). 

0 : Non-dimensional fuel power density. 

Q” : Non-dimensional heat flux, 

h : Heat transfer coefficient (W/inV 

2 ^ 

: Steady state heat transfer coefficient (W/m / C). 

H : Non-dimensional heat transfer coefficient, 

k : Thermal ccnductivity of cool^t (W/m/ *C). 

L : Length of fuel rod (m). 

Re : Reynolds number. 

Pr : Prank nimiber- 
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Coolant flow rate (Kg/s). 

Steady state coolant flow rate (Kg/s). 
Non-dimensional perturbation in coolant flow. 
Cross-section area for coolant flow (in )- 
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coolant density (Kg/m ). 

o 

Specific heat of coolant (J/Kg/ C). 

Perimeter of the fuel rod (m). 

Non-boiling height (m). 

Steady state non-boiling hei^t (m). 

Perturbation in non-boiling height (m). 

Coolant flow velocity in two phase region (m/s). 
Steady state two phase cool^t flow velocity (m/s). 
Perturbation in two phase coolant flow velocity(m/s) 
Coolant flow velocity in single phase region (m/s). 
Steady state single phase coolant flow velocity (m/s) 
Perturbation in single phase coolant flow velocity. 
Characteristic frequency of phase change i/s). 

Chang© in specific volune of coolant at phase 
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change (m /Kg). 

Specific volime of coolant in liquid phase (m /Kg) 
Latent heat of coolant (KJ/Kg). 

Enthalapy of coolant in licMid phase (KJ/Kg). 
quality of coolant. 

Steady state quality of coolant. 

Perturbation in coolant quality. 

Void fraction. 

Perturbation in void fraction. 

Perturbation in core average void fraction. 
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Cxx)lant density in liquid phase (Kg/m ). 
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Coolant density in vapor phase (Kg/tn ). 

3 

Density of two phase mixture (Kg/m ). 

3 

Stscidy state density of two phase mixture (Kg/m ). 
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Perturbation In two phase mixture density (Kg/m ). 
Time required for a fluid particle to loss its 
subcooling (s). 

Enthalpy of coolant at core entrance (KJ/Kg). 

Axial coordinate referenced to the core inlet (m). 
Axial coordinate referenced to the boiling 
boundary (m). 

Time varicbl© (s). 

Time at which fluid pau'ticle crosses the boiling 
boundary (s). 

Time at fluid particle reaches at the top of 
core (s). 

Time vari^le referenced to t^ (s). 

Perturbation in t (s). 

Time recpiired by fluid particle to pass from the 
two phase region (s). 

Friction factor. 

Hydraulic otean diameter (m). 

2 

Gravitational acceleration (m/s ). 

Pressure drop in single phase region (W^). 
Perturbation in pressure drop in single phase 
region (MPa). 

Pressure drop in two phase region (fff'a). 
Perturbation in pressure drop In two phase 
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region (MPa). 

P : Delayed neutron fraction. 

: Decay constatfit of precursor (/s). 

A : Average life of prompt neutron (s). 

$ : Non-dimensional reactivity of fL»l material. 

: Void reactivity coefficient ($/ Z void). 

DR : Decay ratio. 

Ap : Total pressure loss in core (MPa). 

loss 

r Total frictional pressure loss in core (MPa). 

Ap : Total acceleration pressure loss in core (W>a). 

a 

Ap : Total pressure loss due to resistance to flow at 

c, * 

abrupt are changes (MPa). 

Ap : Total driving pressure (W>a). 

u 

3 

: Coolant density in down comer (Kg/m ). 

dc 

3 

: Coolant density at channel exit (Kg/in ). 

X : (i»lant quality at exit. 

ex 

3 

p : Average coolant density in channel (KgA* )• 

Ljj, : Riser height (m). 

NOTE : Terms having bar(”) at top of any variables represents 
the Laplace Transform of that variable. 



CHAPTER 1 
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INTROOUCTION 


The 1990s' has witnessed the spawning of new nuclear reactor 
designs In US and Japan which span the technological range from 
evolutionary light water and heavy water, through light water 
passive and gas cooled, to llcjild metal fast breeder. 

The likelihood of successful development and operation of any 
Indivickial reactor design depends on a mix of both * hard ' and 
'soft ' variables. The hard variables are the more objective 
factors such as actual design characteristics, cost, length of 
technology leap, and degree of research and development support by 
the government and Inckistrlal organisations. 

The soft variables are more ambiguous and depend more on so 
called ' public perception ’ rather tlwi on facts. They may 
Involve Imponderables like attitude towards carbon dioxide 
production and global warming, and alternative power technologies. 
Other factors are more nation specific. Including degree of 
political opposition to siting of nuclear plants or radiation 
waste storage sites, public ambivalence tousards the unknowns of 
the new technology versus the fami liar one and the regulatory 
proceckires. Yet other soft variables are utility specific such as 
size of available markets and suit^Hlty of a design's power 
output - 

Most of the designs developed so far have aimed at Increased 
safety and economical operation by sln^jllfylng construction, 
operation and ma1ntenax« because these designs have Included 
these factors as principal objectives. 
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The Advanced Bolling Water Reactor (ABWR) and its smaller 
counterfiart the Simplified Bolling Water Reactor (S8WR) are among 
the newer designs In BWR proposed by General Electric (GE) 
(Reference New Reactor [6])- The ABWR and SBWR designs are expected 
to show Improvement over the current fleet of BWRs In plant 
availability, operating capacity factor, safety and reliability 
while reckiclng power generation costs, construction times, 
occupational radiation exposure and radioactive waste. By 
simplifying desl^ of the components system and structure »id by 
using natural circulation of coolant further Improvements In 
safety, performance and economy can be made. 

The designs of ABWR (fig 1.1) £tf>d SBWR (fig 1.2) have much 1n 
comnon, the only difference between them are the power rating, 
core flow recirculation (10 Internal punps for ABWR and natural 
circulation) and extent to which some of safety systems uses 
active versus passive features. An In^Kwtant chau^terlstlc of 
both ABWR and SBWR Is the elimination of external circulation 
piping and therefore permits a compact containment desipi. It also 
allows elimination of large vessel nozzles below the core, and 
therefore the desl^ of a more economic emergency core cooling 
systems. Elimination of external circulation piping results In a 
greater then 50 X reduction In welds, less In service Inspection 
(ISI) of the primary system boundary, less occupational dose 
during ISI, increased system integrity and sizable saving 1n 
C 2 ¥> 1 tal cost. 

In response to the Increasing Interest In potential future 
nuclear units comprising the characteristics of smaller size, 
greater simplicity snd more passive safety features, GE started 
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studies In 198? of a 600 MW© BWR with simplified power generation, 
safety and heat removal systems. The basic objectives that were 
established for this new design SBWJ are 

► Power generation cost to be cheaper than those of coal. 

► Plant safety design simpler than 1n current design by using 

► passive safety concepts, 

► Design based on existing technology. 

► Shorter construction schedules. 

► Electrical rating 1n the 600 MWe range. 

► Improve the seismic resistance of the core. 

► Length®^ the continuous operating period. 

► Simplify the system by eliminating movable components. 

► Improve the operability and the maintainability. 

Selection of the natural circulation as the means for 
providing coolaff^t flow through the reactor, coupled with a 42 Kw/1 
core power density, results In a number of benefits to help 
satisfying SBWR objectives, Qx^jared to the existing, forced 
circulation plants, the natural circulation BWR offers low fuel 
cycle costs, fewer operational transients, and Increased thermal 
margin for the transients expected to occur. In addition, 
elimination of the recirculation, pumps and controls needed for 
forced circulation substantially simplified the design, (fig 1.3) 
Conventional BWR safety / relief valves - which opens and 
closes to discharge reactor vessel steam to suppression pool - are 
not needed In SBWR because In It an isolation condenser Is placed 
In the Isolation condenser pool so when the reactor vessel Is 
isolated from the turbine condenser, the Isolation condenser 
controls the reactor pressure automatically without the need to 
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remove fluid from Uie reactx>r vessel. 

Gravity-dr Iven core cooling system provide a simple approach 
to emergency core cooling system because It eliminates the need of 
pumps. Simplification of the components ^d the system made the 
primary containment vessel and reactor building smaller which 
reduces the radiation exposure and construction cost 
substantially. 
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CHAPTER 2 

PRINCIPLE OF NUaEAR REACTORS 


2.1 Introduction: 

The thermal energy procfciced in a fission power plant is the 
kinetic energy of the fission fragments and to a lesser extent of 
the emitted neutron and other particles and radiation such as 
gamma rays which get converted to heat when these particles are 
absorbed. This heat is removed by a coolefit and subsequently 
utilized in a thermodynamic cycle. 

hkjtclear reactors are variously classified according to 
general purpose or function, type of the moderator, type of the 
coolant, neutron energy classification, type of fuel, type of 
core internal design and other factors. 

While the other possible coolant besides light water can be 
used like the heavy wetter, organic liquid, gases and liquid 
metals, light w^ater has been used most extensively , becaise of 
its availability, low punping pcwer and the advaxed state of 
knowledge concerning its chemistry and thermodynamics properties 
etc. 

2.2 Lij^^t water Reactors: 

The reactor in which the coolant is light water may be 
broadly divided into PWR and BWR. As the names suggest. Boiling 
water Reactors are those in wrfhich the woter boils within the 
reactor, and in Pressurised Water Reactor the pressure is high to 
avoid boiling. Thus the water pressure corresponds to the 
saturation tenperature at the reactor pressure in a BWR. The 
reactor pressure is roughly between 600 and 1000 psia. A schematic 
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arrangement of a BWR is shown in the f1gure( 2.1 ) 

In a Pressurised Water Reactor, reactor pressures are hl^ier 
of ttie order of 2000 psia, for the primary coolant loop within 
the reactor. A heat exchanger in which the primary coolant 
exchanges heat with a secondary coolant is used. A schematic 
sketch of a PWR is shown in the figure( 2.2 ) 

However, in this study, we focus our attention on the former 
type of reactor that is the BWR. The boiling reactor has a 
function closely resembling that of boiler In a conventional 
fossil-fuel steam power plant and is basically simpler than it. 
While in a boiler heat is transmitted from the furnace to the 
water indirectly - partly by radiation, partly by convection and 
partly by conduction; with oonibustion gases acting as an 
intermediate agent or coolant, in the Boiling Water Reactor, the 
coolant is in direct contact with the heat producing nuclear fi»l 
and boils in the same compartment in which the fuel is located. 

The sinplest form of a boiling reactor power plant as shown 
in the figure( 2.3 )consists of a reactor, a turbine generator, a 
cond^ser and the associated equipments ( such as ejector, cooling 
system etc.) and a feed pimp to force the Incoming coolant to 
core. Slightly sub-cooled liquid enters the reactor core in the 
bottom where it received the sensible heat for saturation plus 
latent heat for vaporisation. When It reaches the top of the core 
it has been converted to a very wet mixture of liquid and vapor. 
The vapor separates from the liquid, flows to the turbine, does 
work and is condensed by the condenser, and the condensate is then 
piBn^ed back to the reactor by the feed pump. That is how the power 
is generated in the conv«^tior>al BWR. Such a cycle is called 
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Schematic arrangement of liquid-cooled-reactor power plant. 
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direct cycle BWR, whidi has the disadvantage of primary coolEff>t 
coming in contact with turbine etc. To overcome this problem one 
can employ Indirect cycle BWR In which a heat exchanger 1s used to 
separate the radioactive primary coolant from the non -radioactive 
secondary coolarit. 

2.3 Natural Circulation BWR: 

The simplified and modified version of the conventional BWR 
is the natural circulation BWR, more generally known as Simplified 
Bolling Water Reactor (S8WR). The present thesis deals with this 
simplified BWR. Except for the neutronic / thermal-hydraulic 
des1{y>, there is no difference between the desl^ of natural 
circulation BWR and forced circulation BWR. In the natural 
circulation BWR the necessary recirculation flow fore core cooling 
will be produced by the driving force originating from the 
tftioyancy dute to density head d1fferer>ce In the closed loop of the 
recirculation flow. (In the absence of recirculation pums) 

Actually the technology for the natural circulation Is not 
new to BWR. The Dodf^waard plait in Netherlands has operated on 
IMs principle at a life time capacity factor of 84Z. Larger BWRs 
have be^ operated at 501 power levels in natural circulation mode 
to prove that BWR of this type Is Indeed possible. 

A schematic of the natural circulation BWR as It contrasts 
with the forced circulation BWR Is shown In the fIgureC 2.4 ). 
The major design aspects of the natural circulation BWR has been 
shown in the figure ( 2.5 ). 
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Fig.2.4 Schematic of (a) forced-circulation BWR and (b) 
natural-circulation BWR. 
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Figure (2-5): Major desiffi aspects of natural circulation BWR 
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CHAPTER 3 

PROBLEM STATDCWT AND PROCEDURE 

3. 1 Problem Statement: 

The present study investigate the feasible power rates for a 
natural circulation BWR. For this the core stability of the 
reactor was to be examined. It further studies effects of the 
variation in active fuel length and core power density on the 
necessary riser height for a natural circulation BWR. 

3.2 Procedure : 

In the analysis, a simple flow model was adopted where the 
main recirculation flow in the riser, downcomer and parallel core 
channels was assLsned to be vertically one dimensional. 

A code was developed to analyze the nuclear-coupled stability 
as well as thermal-hydraulic stability of the fuel channels. This 
st^ility code consists of a reactor point kinetics model and an 
axially one dimensional thermal hydraulic model of the fuel 
channels, and it adopts approximately the same line of approach as 
used in L^»y & Moody [?]. This stability code calculates the 
frequw>cy response of the open loop of the system and obtains the 
decay ratio. The decay ratio is defined In chapter(4) equation 
(123). The input parameters for this code are n&itron data as well 
as the thermal hydraulic data and the channel flew. 

Finally by balancing the total pressure drops in the 
respective parts of the recirculation flow loop formed by the down 
comer-core- riser with the total driving pressure, the necessary 
riser height is calculated. 

figure( 3.1 ) shows the procedure for calculating the 



17 


necessary riser hel^t. First we obtain the minimum recirculation 
flow rate necessary to satisfy a tentative requirement on the 
decay-ratio (DR) as 

DR < 0.25 

This requirement for decay ratio Is conservative enough to 
give an overall Inherently safe reactor design. Finally based on 
this mlnimim recirculation flow rate necessary riser hel^t Is 
calculated. 



FlgureC 3.1 ): Calculation procedure for iwcessary riser height. 
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THEORETICAL ANALYSIS 


A.l Types of IristabUUles In BWR: 

4-1.1 Nuclear coupled Instability: 

During the early days of BWR technology, there was 
considerable concern about nuclear-coupled Instability ; that is 
interaction between the random boiling process and void- reactivity 
feecftDack modes. Indeed, it has been proved by extensive series of 
experiments that while instability was observed at lower pressures 
but it was not expected to be a problem at the higher system 
pressures. In natural circulation BWR there is some probability of 
arriving at low reactor pressure hence nuclear coupled instability 
affialysis become a necessary for this reactor. 

In addition to nuclear coupled instability, there are number 
of static and dynamic Instabilities. Some of the static and 
dynamic instabilities which are generally considered in design of 
BWR, are 

4.1.2 Flow Excursion Instability : 

It is concerned with the interaction between the pu^'s 
head-flow characteristics and the hydraulics characteristics of 
the boiling channel. Since in natural circulation pumps are 
eliminated hence this instability does not occur. 

4.1.3 Flow Regime Relaxation Instability : 

It is due to the change in the flow regimes. This instability 
is self stable in case of natural circulation BWR because if due 
to a small perturbation a slug flow pattern changes to asinular 
flow a lower pressure drop is experienced which may tend to 
increase the flow rate and thus cause the flow regime to return to 
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Hs previous state. 

4.1.4 Density Wave Oscillations : 

This Instability Is due to the feetfcack and Interaction of 
various pressure drop components and is caused specially by the 
lag introcixced throu^ the density head term due to the finite 
speed of proportion of kinematic density wave. This instability 
is obviously quite important and needs to be examined. 

4.1.5 Pressure Drop Oscillation : 

It is the instability mode in which flow excursion 
Instability and a compressible volume in the boiling system 
interact to produce a fairly low frequency oscillations. This 
Instability is normally not a problem in natural circulation BWR. 

There au'e however, other instability modes that can occur in 
two phase systems. Since here we are concerned exclusively with 
natural circulation BWR technology, we can neglect them here. 

4.2 Stability Analysis: 

4.2.1 Principle: 

Consider a subcooled flow entering a heated channel as shown 
in figure ( 4.1 ). The boiling boundary, X(t) , defined as the 
instantaneous location of the point where the bulk fluid 
temperature reaches saturation temperature, oscillates due to the 
inlet flow and power density fluctuations. Change in flow and the 
length of the single phase region conbine to create an oscillatory 
single phase pressure drop. At the boiling boundary these 
perturbations are transformed into quality (or void fraction ) 
perturbations that travel up the heated channel with the flow. 
Combined effects of flow and void fraction perturbations and the 
variation in the two phase length create a two phase pressure drop 
perturbation. However the total pressure drop across the boiling 
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clwtnel »-»as to remain constant with time in BWR, thus, the two 
phase pressure drop perturbation produces a feedback perturbation 
of the apposite sign in the single phase region. 

To simplify the analysis and yet retain the essential 
features, the following assumptions are considered in the 
present stability amlysis of natural circulation BWR. 

• Diabatic homogeneous two phase flow. 

• Constant system pressure. 

• Uniform axial heat flux. 

• No subcooled boiling. 

The basic solution scheme consists of coupling analytical 
models for the single and two phase region with a model for the 
dynamic thermal response of the heated surface, such that imposed 
boundary conditions are satisfied. 

4.2.2 Model For The Fuel-rod: 

The present analysis for fuel rod essentially follows the 
treatise given in Lahey & Moody [73. 

Using a single-node luiped heat capacity model one cat) write- 


M-C = q-V - q"*A 

' dt 


( 1 ) 


where 


M 

C 

T 

q 

q 

V 

A 


f 

f 


II 


= Mass of the fuel rod. 

= Specific heat of the fuel rod. 


= Fuel temperature. 

= Fuel power density. 

= Heat flux. 

= Volume of fuel rod. 

= Surf^ace area of fuel rod. 
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At steady state 



q^-v = q;;-A 


( 2 ) 
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where 


= Steady state power density, 
q'^ = Steady state heat flux. 

Let 



Q" 



(4) 


where 



(5) 

( 6 ) 


Steady state fuel teoiperature. 

T = Coolant temperature. 

C 

Steady state coolait temperature. 

Eliminating q, q" and from (1) by using (3), (4) and (5), 


M*C 



Q.y.q^- Q".A*q“ 


(7) 


Now taking Laf>lace Transformation of ^x>ve equation, 

M-C -s-Ii^ = q„*V*® - q"*A*Q" (8) 

f f 0 0 

Now to express perturbation In heated wall temperature 
terms of heat flux perturbation ( Q") use - 


q" = h-(T -T ) 

1 C 

where 


(9) 

( 10 ) 


in 


h = Heat transfer coefficient of coolant. 
h^= Steady state heat transfer coefficient of coolant. 
From (4), (9) and (10) and after simplifying. 
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^ - ho h 

Q = + ( n - n ) (11) 

h a" 

0 ^0 

The second term in ri^t hand side in the above ecpjation is 
non-linear it can be linearized by replacing h by h^ in that term 
hence 

h 

Q" = H + -^*(n - n ) (12) 

q.. r c 

here 



Now from Reynold's Analogy 

K - n 03 ^ rv-.0.8 _ 0.3 

h - 0.23’ ’Pr 

where 


(13) 


(14) 


k = Thermal conductivity of coolant. 
L = Length of fuel rod. 

Re = Iteynold's Nurter 

I 

= constant*G 


G = Coolant flow rate 
Pr = Prandtl Number 


Since k,L and Pr are constant for given cool^t and operating 


pressure of system hence 

h = constant ’G^ ' ® (15) 

and h = constant*G°' ® (16) 

0 0 

where G = Steady state coolant flow rate. 

0 

Differentiating (15) with respect to G, 

constant*0,8’G (17) 

clG 

Combining (16) and (17) , 
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cfri 


cfe-G 


- 0.2 


= 0 . 8 ' 


. 0.8 


(18) 


To linearize equation (18) replace G by , 


dh dG 

— * 0.8 (19) 

h, G 

0 0 

or H = 0.8-C (20) 

here € = non-dimensional perturbation in coolant flow rate 
Substituting H from (20) into (12) and taking laolace 
Transformation of the resulting equation one gets 


Q" = 0.8-C + — -(D - n ) (21) 

n f c 

^0 

Now eliminating from ( 8 ) and (21) and sln^jllfying using (2) it 
reduces to - 


M-C -s . r 1 

Q.. s f L ]. 0.8-C + — - — C - — - R — (22) 

M-C^-s + h„-A *- M-C^-s q" ® 

f 0 f 0 

Equation (22) gives the response of the heated wall in the single 
phase region to perturbations in coolant flow, core power daisity 
and coolant temperature. 

4.2.3 Model for The Single Phase Region: 

The analysis for the single phase regicn essaitially follows 
same line of approach as used in Lahey & Moody [71. 

Single phase energy balance can be written as 
57 9T 

^ ,p .0 . — ? + G-C • — ^ = q"-A (23) 

c c c ^ c « p 


where 

A = Cross-section area for coolant flow. 

C 

p ' Coolant density. 

C 

C = Specific heat of coolait. 

C 

A = Perimeter of fuel rod. 
p 
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t arid 2 are time and distance (from core Inlet) varidales. 
At the steady state equat1on(23) reduces to. 



q" - A 


— 

1 = ° P - 

° G -C 

0 c 

(24) 

^ dz ^ 


Now substituting 

T = n + T „ In 

c c cO 

(23) and then simplifying , 

an 

an 

r ^ ^ 

A -p -C • — - + 

G-C + G-C 

- — -] = q"-A (25) 

c '^c c at 

c ^ c 

oz 

'^az ^ 


M^ing the above equation linearize and subtracting from 

both sides of equation (25) , 


an an r ^ > 

^ ]o- PSV - -’S’-A. 


Sinpllfying (26) using (24) and transforming q" and G by Q" and 

C , 

A -p an an q"-A^ 

. £. + = _i_L. . (0 - C) (27) 


an 

C . 


' + 

at 



Now talcing Laplace Transform of above equation, 

A ‘p cjn d’'*A 

_£_c .g.n + — S_ = _1_P_ . (0 - C) 

G_ dz G„-C^ 


. (0 - C) 


dZ 

0 C 


Now substituting Qi" from (22) In above ecjiatlon. 


dn _ _ 

X(s)-n + — - = Y(s)*C! + Z(s)-€ 

dz 


where 


r ^ 1 ____ 

x(s) = ‘S + * 1 

^ r G -C M-C -s + h *A J 

'^0 0 c f 0 

,».A r 0.8-M-C -s 1 

0 p . _i 1 ' 

J -C ^ M-C -S + h *A ■* 

0 c i ^ 

1 (32) 


Y(s) = 


Z(s) = 


0 P 
0 c 
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Ec^iation (29) is linear differential equation of first order hence 
can be solved easily, solving it with the assun^tion that 
perturbation in inlet coolant ten^rature is zero, one gets 

1 - 

n = [ Y(s)-G' + Z(s)-« l-z- [ 1 (33) 

^ *- X(s)-z -1 

Next step in the analysis is to determine the dynamics of the 
boiling boundary- For the assumption of uniform axial heat flux, 
the coolant enthalpy increases lineau'ly with axial distance from 
the bottom of fuel rod- In fact a positive perturbation in cool^t 
temperature at the boiling boundary causes a negative perturbation 
in the boiling boundary as shown in fig. ( 4.2 ) 

Let X = Non-boiling height 

X^ = Steady state rton-boiling height 
X' = Perturbation in non-boiling height 
Now making energy balance in the length X' at boiling boundary 
X' -A -q" * - G„*C -n ( X^) (34) 

p 0 0 c c 0 

Taking Laplace Transform of above equation , 


!• = — . n ( X ) 

Q"-A ' ^ 

0 p 


(35) 


Now substituting R ( X ) from (33) into (35) one gets 

c 0 


G *C r 

X' = - — ? • [ Y(s)‘C + Z(s)-0]-X^* 

q" -A *■ 

0 p 


1 - exp( X(s)*X^) 


X(s)*X, 


( 36 ) 


This can be written as 
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h 


h 


f 



Fig. 4.2 Boiling boundary perturbations. 
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X' = A^(s)-^- 

+ A^Cs)-® 

(37) 

where 




1 - exp( X(s)*X ) 


A (s)= - ^ . 

0 

• Y(s) (38) 

q" 'A 

0 P 

X(s)-X^ J 


and 



G -C 

A (s)= - . 

■ 1 - exp( X(s)-X^) n 

■ Z(s) (39) 


K(s)-X„ J 



The main purpose of the single phase analysis has been to 
derive the expression for the dynamics of the boiling boundary for 
which equation (37) is the required result. 

4.2.4 Model for Two-Phase Region: 

The analysis for two phase region is similar to as done in 
Lahey & Moody [7]. 

For constant system pressure, the continuity equation for a 
homogeneous two phase mixture can be written in the form 
a"* A *7 

= __0_P__£3_ ^ 

^ h -A 

fg X 

The parameter , Cl .which has the units of reciprocal of time, 
is frequently referred as the chairacteristic frequency of phase 
change. It physically repres«it the speed at which phase change 
takes place - 
here 


V = coolant flow s^elocity. 

V = change in specific volume of coolant at phase (^lange. 

I'g 

h = enthalpy change of coolant at phase chansp. 
f 9 

A = cross section area of coolant flow. 

X 

Corresponding energy equation for a homogeneous two phase system 


is given by 
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DX A r. 

_ - fi-X = 0 

here 

X = qualHy of coolant 

v^= specific voluTte of coolant In liquid phase. 

Equation (40) can be integrated along constant time characteristic 
to yield, 

Z 

v( 2 ,t) = V (t) + / ft'dz' (42) 

X(t) 

here 

v^(t) = single phase inlet flow velocity at time t. 

X(t) is the instantar^eous position of the boiling boundary which 
can be written as 

t 

X(t) = / v^(t')*df (43) 

t-u ' 


The parameter p is the time rajiired for a fluid particle to lose 
its sub-cooling. It can be evaluated by making energy baleaxe in 
the sub-cooled section, i.e. 

q''*A ‘L •t' = (p ’A "L )*(h - h ) (44) 

p* cxs fi 

here 

L = Length of sub-cooled section. 

8 

h = Enthalpy of coolant at core entraaice. 

i 

hence 
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= v^(t) 4 f ft.dz- ( 47 ) 

X(t) 

W " -* ~ (48) 

^ ~ (49) 

The next task Is to Integrate eq" (49) to obtain the locus of 
the fluid particles In space-time plane. It Is convenient to 
transform from the axial coordinate, 2 , referred to the core 
Inlet , to y, which Is referred to the boiling boundary. 

Thus y = 2 - X(t) (50) 

Differentiate (50) and (43) with respect to t , 

^ - dX(t) 
dt dt dt 

and 

= v^(t) - v^(t - P) (52) 

Now combining (49),(50),(51 ) and (52) one obtains 
- 0-y = v^(t - P) (53) 

This Is a linear differential equation of first order hence caff> be 
solved as 


t 

y(t) = exp(fi*t)* / exp(-ft‘t' )‘Vj(t'- p)'dt' (54) 

Here t^ Is the time at which fluid particle crosses the boiling 
boundary. 

Let 


t = t - t^ 


t = 

1 



(55) 


dr^= dt* j 

using (55) In (54) one gets 

T 

y(r) = exp(fi*T)- / exp(-ft-t:^)-v^('r^- t^)*dr^ (56) 

0 

Equation (56) can be perturbated with respect to v and r and 


Laplace Transfonned to yield 
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and 

here 


^ _ 

dx 


exp(-s-i^)-[ 1 - exp{(fi - s)*t}] 
(s - Q) 

exp(O-t) (58) 


(S7) 


Steady state single phase coolant flow velocity. 

Let 


T' = perturbation in r 
and = perturbation in 

Now T' must be expressed in terms of C and C , for this 
perturbate and Laplace transform equation (50) at a fixed axial 
position ,2 , then 

^ ay 

dz = 0 V’ + t' + X' (59) 

av^ ar 

By combining (57), (58) and (59) one obtains 

[ X' + .( 1 - exp{(Q - s).T} )• v; J - (60) 


But 


so 


G = V -A -p 

1 X c 


c = 


dG 


dv , 


10 


hence v’^ = C 


(61) 


Now substituting from (61) into (60) , 

— i) • T r S • ^ > 

[ ^' ^ fi^ ’ ( ^ J ’ '' 


r*= - 


10 


(62) 


To relate perturbation in coolant cjiality with perturbation 


in t one proceeds as follows- 

Solving equation (41). which is a linear differential equation of 
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>:(t) = e 


fi-t 




(63) 


Using (5S) the equation (63) can be integrated and simplified as 

V 


x(x) = ^ .( - 1) 

f g 


(64) 


Now by perturbating and Laplace Transforming equation (64) one 
gets 

V 


X’ 


f 0-T 

e • T' 


(65) 


fg 


here x* is the perturbation in the coolant cpjality. 

In the BWR, the nuclear feedback is largely through void 
reactivity coupling, thus ,it is important to relate quality 
perturbations to perturbation in void fraction. A reasonably 
accurate void-quality model has given the relation between void 


fraction and the quality as- 

X 


a 


r Pa ^ 

X + *(1 - 

L Pf J 


( 66 ) 


where 


a = void fraction 


p »= coolant density in vapor phase 
g 


p = coolant density in liquid phase 
f 


and Cq = t (0.71 + 0.0001 -p) - U 


-1 


here p = pressure of reactor (in psia) 


By perturbating and Laplace Transforming equation (66) one gets 
«• = K * X’ <6"^) 

where 


K = 


C 'P • 
0 


i2 


( 68 ) 


P 
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here 

a' ' perturbation In void fractiori 
x^= steady state coolant quality 
By combining (62), (65) and (68) one gets 

1 

«•= - K-0~- — + - exp{((^s)‘T}]-€ — (69) 

g *- 10 J 

Now to get perturbation in the core-average void-fraction defined 

as 


a* 

• V g 


J a' -dz 


-j-* / a' - dy 
0 


f- I 


(t -t^) 

ex 0 


«• 


dy 

dt 


0 

X 


-1 -L 

1. r 

L J dt 


dr 


0 

r 


i- f-- 


'^10* ® 


ft-T 


dt 


dT 


(71) 


0 

Here t is the time required by fluid particle to reach at the 

ex 

top of the heated wall, and r is the time in which the fluid 
particle has travelled in the two phase region. 

Substituting a- from (70) into (71) and solving the resultant 


equation. 


a' 


K*v, 


L*v 




X' • [ exp((l‘t ) - 1] + B(s) 

ex 


(72) 


where 


K'O'V 




— V 


r 


exp(ft*t )-l exp[(2-(l-s)*t^^]-l 


10 


(s - 0) 


L' 


B(s) » - 


t « u 


0 


(2*0 - s) 
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Making first order approximation to the exponential function in 
expression of B(s) theri we find that It will become zero. 
Her»ce after making first order approximation In (72), 


K-v, 


a' 

avg 


Q • r 


L 'V 


(74) 




Now substituting for X' from equation (37) 1n (74), 


a' 


avg 


T^(s) • C + T (s) - « 


(75) 


where 


T^(s) = - A 


,(3). 

L*v 

fg 


(76) 


and 

r K*V ‘Q'T 1 

l^{s) = - AgCs)- ^ (77) 

■- L*V^ -* 
fg 

Now substituting Aj(s) and from (38) and (39) Into (77) 

and making first order approximation to the exponential 
function. 


6 -C 

T^(s) = Y(s) 


q" -A 

0 p 


r K*V -1 

L L.y J 

fg 


(78) 


and 


6^-C 

^0* Z(s) 


q*' -A 

0 p 


K*v 


•it' L -I 

f ex 

L-v. 


(79) 


fg 


Now substituting Y(s) and Z(s) frcm (31) and (32), 


r K'V,*0*t 1 r 0,8*M*C *5 -j 

\ ^ 1 ■ • ■ ' J 

L L-v. J “ l-H.c.sth-A J 


T^(s) » 


(80) 


fg 


and 


r K'V -O-t 

f e’* 

1 • X . [ 1 

(81) 

fg 

■1 ° L M-C^'S + h^-A •* 
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X = 


J_ 

Q 


In 


[ ] 


( 82 ) 


henoe 


X = 
e X 


i- ^[ i)-(L-V , v„ ] 

0 10 ■* 


0 ' 10 

Equation (75) gives the perturbation In the void-fraction tije to 
perturbation in coolant flow and core power density. 

Now to analyse the Instability due to density wave 
oscillation, perturbation In the pressure drop In single phase 
and two phase region are to be known. 

4.2.5 Analysis for single phase pressure drop perturbation : 

The analysis for this section essentially follows treatise 
given In Lahey & Moody [7]. 

The appropriate momentum equation for the single phase 
portion of the heater can be written as 


' ■■■ = D • 


dv, 

1 cl 


dz 


at 


2-D, 


+ p^. g 


(84) 


where 


f = friction factor 

D = Hydraulic mean diameter of coolait flow 
h 

g = Gravitational acceleration 

In terms of the coolant flow rate equation (84) can be written as 


JP. - _L . + 

■ az ■ A dt 


f-G" 


2.D^*A^-Pc 


_ + p^* g 


(85) 


Integrate (85) with respect to ,z, frcm ,z = o, to ,z = X, to get 
the pressure drop In single phase ,(p^)» QS 


p - / f ■ -f- f ® ]■* 

0 L 2*D 'A *P 

Tl» perturbation In slngla phase pressure drop .(P,). can be found 
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by differentiating (86) with respect to G and X and then adding 
both. After linearizing the resulting equation It reduces to 

p . 1 ^ , r ^‘^0 ^ 

1 J L A dt ^ * 5 — + P • g -x* 

0 -A .p -I L 2-D -A^-p ^ J 

^ ^ h X c 


3 -A -p 
n X c 


Replacing 6 g by G^- C and then taking Lapla» transforming In 


above equation. 


'f’G^’X T r f‘G~ 

-•s-G^ + 0 J. c + ^ 0^ 

X D -A *p -■ L 2*D -A "p 

“ ^ c h X c 

Substituting X* In terms of C C from (37), 


+ p^. g 


= r^(s) . (E + T^{s) • il 


where 


r/s) = 


r_^ 

^ 2 *D • A^* 


1 X^-G^ 

— r — + p • g ]• A (s) + — + — r — 's 

*A *P D "A *0 X 

h ''x ''x * — ( 90 ) 


r3(s) 


r 1 

= z — + 9 

L 2-D -A^-p ® J 2 

h X c 


(91) 


Hertce resultant total single phase pressure dr<x) perturbations 
transfer functions can be written by (89). 

4.2.6 Analysis for two phase pressure drop perturbation : 

The analysis for two phase pressure drop perturbation 
essentially follows approach as given 1n Lahey & Moody [7]. 

For the case of vertical, homogeneous two phase flow, the 


momentum equation is given by 


dz 



(92) 


Where p 1s the density of two phase mixture which Is given by 
H 
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p = 

H V,+ X*V 

f fg 


SubstUuting x from (64) 1ntx> (93) simplifying. 


Ph = Pf-e (94) 

lo get two phase pressure drop (p^) Integrate ecpiatlon (92) 
with respect to z from ,2 = X, to ,z = L . Hence 


^2 ' I h' 




* P„- B • * 
. 


Perturbation In the two phase pressure drop (P^) can be found 
by differentiating p^ with respect to v and X thai adding. 
Finally taking the La^jlace transfer of resultaa^t expression. 


. r_l 5 

^ 2-D. -j 


4 a]- 5' 

^ -A “* 

h X c 


The next step 1n the analysis is to specify the various terms 
In the Integration of (95) so that spatial Integration can be 
performed. 

First consider the two phase velocity perturbation, (v') , - 
Equations (46) and (49) can be combine to give 


v( 2 ,t) = v^(t) + 0*[z - X(t)] 


After perturbating at a fixed z and Laplace transforming. 


= - ft-X’ + V' 


0-X' + 


p ’A 

C X 


Next term to be considered 1s the two phase acceleration 

perturbation , ^ j , - 

combining equations (50) and (97), 


j(yj t'l = V ^t) + Q* 
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Now taKUig material derivative of both sides of equation (99) and 
conil)1n1ng the result wlUi equation (53), 

dv 

= 0-[0-{2 - X(t)3 + v (t - y)] + (100) 

Dt 1 dt 

Perturbation In [-5^] at fixed 2 will be given by 


= _ o^.x' + 0-v'(t - y) + -jr^ (101) 

I Dt J 1 dt 

Now taking Laplace transform 

r.25L.y = _ 0^-X' + D - + s-v' (102) 

I Dt 3 1 ^ 


using (61) 1n (102) 




+ f Q 


a. ^ 

e + s 


j p •, 


p ‘A 

c X 


(103) 


Tbe next term that must be considered is the two phase mixture 
density perturbation , (Pjj)- 

Perturbating equation (94) and Laplace transforming, 
i' = - O-p (104) 

using (62) in (104) 

2-ft-T 


ft-p^*e 

-¥ 

^ 


L f.5 < , c. 


Integrating (58) and combining it with (99). 

Q-v 

v^(2) = {>*(Z - V ^ ''to ' ''10* ® 


( 106 ) 


So 


= '"to* 


(107) 


Combining (94), (96) and (106) yield 


''to-^f 


HO + 


ft-(z - X„) 


( 108 ) 
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Nc^ using (108), (107). (105), (103) and (98) equation (96) can be 
Integrated then one obtains 


^2 = n/s) • C + 11 (s) . C 


(109) 




n, 


(S) . A (s)-[ ( - n-t ) - .t,.(L - X )1 


and 


( 110 ) 


II 


,(s)-A,(s)-[(- ^°° 

Z’D ‘A ‘P 


f -G. 




h-x ^ 


A -D, 

X h 


-v] 


f-p 


+ p -T -O-Cl - S*IA) + p ‘t *5 + 

C ex ex D 


•Q-(L - X^) 


o n c.i)\ *p^*g*t 

+ 2-p -0^. .11 : 1 ?L_£ 

^ (s-0)^ G„ 


(i- 


(Ill) 


The above equations are adequate to appraise the 
hydrodynamics steAjIllty margins of a boiling system. For 
appraising the likelihood of density-wave-osclllatlon In the 
boiling loop, it Is considered that static pressure is continuous 
around the closed loop, thus 


P + P = 0.0 (112) 

1 2 


hence 


or 


r^. « + r^- « + n,* ® ' 0.0 



^2 -1 

r, - n J 




(114) 


(113) 


This gives the stability criteria for the density wove 
oscillation. 

4.2.7 Neutronic Transfer function (MT) ■- 

A point-reactor model with single delayed neutron group can 


be written as 



4 ) = Nrr(s) • $ 


(115) 


where NT(s) Is given by 
b* (s + 1 ) 


here 


NT = 


b = 


s*(s + b +1) 

P 


(116) 


T^T 


(117) 


and P = Delayed neutron fraction 


X = Decay constant of precursor 

p 


A * Average life of prompt-neutron 
$ = Laplace transformed perturbation In reactivity 


Also, 

$ = K • a’ (118) 

V avg 

where K is the void reactivity feecftiack. 

V 

<h.2.8 Analysis for decay ratio : 

Now fron equations (115), (118), (114) and (74) we construct 


a control loop as- 



The above control loop can be simplify as 
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K -Krr-T 

V 1 

1 - K -NT-T 

V 2 


’’e* \ 

Hertc<» gain transfer function G(s) is given by 


G(s) = 


K -NT(s)*T (s) 

V 1 

1 - K -m-(s)-T,(s) 

V 2 


(119) 


and feedback transfer fur»ctioo H(s) is given by 

rjs) + J\A^) 


Tys) 4 n^ts) 

The overall trarjsfer function will be 

Tfc) T Gts) 

1 - G(s)-H(s) 


( 120 ) 


( 121 ) 


Now applying root- locus technique for stability analysis of 
the system, the poles of the polynomial G(s)*H(s) was found.From 
these ptiles two complex conjugate poles which are closest to the 
origin are most sertsible for stability analysis, were considered. 
Let tl>ese poles be (m + in) 

Tlie behavior of the output with time corresponding to these 
poles will be of like as shown in figureC 4.3 ). 
where 1 is the time period of the oscillation 
Now Decay Ratio is defined as 


(122) 


But T . ^ 
Herrce 
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0 



janod j auo I sustuj p -uou 


eo 

* 


9 


BEHRVIOR OF OUTPUT POWER luith TIME 



^3 


w. - J ( 


173 ) 


^.3 Aivalysls. for rx^ot-ssary riser helc^t : 

irx- afialysis fot rxy';essat y riser heigfil essentially follou's 
lre.<t1es given In El-W.iKil [81. 

Ttie riser Is an untteated extension of the core but usually 
with fewer walls, dividers etc. and of course, with no fuel 
elenrents so that there Is less friction than In the core. The 
basic principle for the calculation of the riser height. In 
natural circulation BWR, Is to balance the total pressure losses 
1n the reel roulatlon loop with the net driving pressure generated 
by the density difference between the water 1n the down comer and 
tliat In the channel. 

The total pressure loss In the core Is given by 


AP 


less 


j: AP 4 E + 1 , 

I a ' 


iUk) 


where 


AP = Total pressure losses. 

loss 

AP = Sum of the frictional pressure losses In the core, 
f 

riser and down comer all computed In the direction 
of the flow 

AP = Sun of the acceleration pressure losses , for this 

B 

000 cxxis 1 cterGcl only in h03t©d chonn©! and 

ignore in the riser and downexjmer since no large 
changes 1n density are encounter there- 
£P = Sum of the pressure losses due to resistances to 

c, e 

flow at abrupt area changes, for this one 
considered only at inlet and exit of the heated 
channel, and neglected the pressure losses due to 
the drag of submerged bodies such as spacers. 


support plates , separators etc. 


Now one wants to know the driving pressure. In case of no 
riser, the driving pressure is given by 


jtaP ~ i hydrostatic pressure in down comer) 
d 

- ( hydrostatic pressure in channel) 

= ( p - p )*L*Q (125) 

d c 

where 

p = density of coolant in the downcomer 

d c 

and p is the average density in the channel and it is given by 
following expression 

n.ico •L +p*L — (126) 

^ ^nb nb b b 

here 


p = Density of coolant in the non-boiling section of the 
''^nb 

channel. 

L - Non-boiling height of the channel. 

nb 

p = Average density of coolant in the the boiling section 
^ b 

of the chanr>el. 

L = Boiling height of the channel. 

b 

p aid p ^ were assumed as equal to (tensity of coolant 
1n the^llquid phase, the eypresslon for \ for the sinusoidal heat 

generation is taken as given- 

NOW to find the effect of riser on the driving pressure it is 
assumed that quality of the coolant in the riser is substantially 
the same as that at the channel exit. If the slip-ratio in the 
riser is assumed to be same as that of in the core, the void 

fraction in the riser will be constant and equal to that at the 

ct^^l exit 1^ the d^sity al^g tl« riser will therefore be 

.... .h... exit. It is then obvious 



45 


tfwil UiP odd'it iofi of a rfser of height .(Lj^) , increases the 
driving pressure by the quaritity (p^ - p )-L„*g . Where p is 

ucoxk ex 

the deiisity of «x)lant at the exit of the channel. 

Thu 5 . the driving pressure in the core of a channel with riser 


is given by 

= t P. •( L + L ) - ( P * L 4 p • L ) ]• g (127) 

G o c K e X K 

The exit coolant density ,(p ), is given by 


= f V, + X -v. 1’^ 

V f ex fg J 


(128) 


and can fje calculated by knowing ,x^^, the exit coolant c^iality. 
Now to get the riser height we equate the driving pressure with 


ttie pressure losses as 

£P = liP (129) 

d 1 08 8 

By substituting from (127) and from (124) one gets a 

equation with riser height as unknown parameter, so it can be 

solved to get the riser height. 


CHW>7ER 5 
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RE SULTS AND DISCUSSION 

5- 1 Var^aHon of Decay RatEo With Core Flow: 

It was found that decay ratio Is generally Improved as the 
recirculation flow Increases, this means that decay ratio 
decreases with the recirculation flow, and the margin of core 
stability Increases. Flgure( 5-1 ) shows an example of the 
relation between the recirculation flow and the decay ratio for 
the two different core power density of 42 KW/1 and 50 KW/l.The 
mlnln«jini recirculation flow necessary for the decay ratio less then 
0-75 Is also shown 1n the same figure. It was found that It 
Increases with the Increases In the core power density. 

5.2 Variation of Minimum Recirculation Flow With Active Fuel 
Length: 

The mlnlmiffn recirculation flew Is found by the condition of 
U>e decay ratio < 0.25 . It was found that It Increases with the 
active fuel length- Flgure( 5.2 ) shows an example of the relation 
between mlnimuBn recirculation flow and the active fuel length for 
tbe core power densities of 42 KW/1 and 50 KW/1. 

5.3 Variation of Riser Height With Active Fuel Length; 

The riser height necessary to produce the mlniraun 
recirculation floe calculated awarding to the procedure set 
out in the Chapter (O. It «s found that riser hel^t Increases 
as the active fuel length Increases. Flguret 5.3 ) shews an 
exaiple of the relation between the riser hel*t and active fuel 
length for core power densities of 42 KW/1 and 50 KW/1. 
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5.4 Variation of Riser Helc^t With Core Power Density: 

As both U»e minimum reclrcxilatlon flow as well total power 

Incrt^ases with core power density hence the 
y r 1s«r helg^jt to provide the minimum recirculation flow 
should Inoroase ai^d Indeed It was found In this study. F1gure(5.4) 
show** Axainple of the relation between the riser height and the 
core powM»i density with keeping cctive fuel length as an 
IndepcHidHnt variable. 

5.5 Conclusion : 

In t-hls thesis the feasibility of the natural circulation BWR 
for higher power rates was studied. The thermal hydraulic 
requlrefnertts for Increasing power rates were Investigated 
analy t Icial ly . It hras concluded that natural circulation BWRs of 
powers > 1000 MW ( electric ) will be feaslbile as far as core 

stability condition Is concerns. The main coficluslons are 
summarised below- 

1 As core flow Increases , the core stability Is Improved- 
► The minimum core flow to secure core stability can be mantalned 
by Increasing the riser height. 

^ The r>©c:essary riser height Is approximately expressed as a 
simple function of power density and the active fuel length. 

In Tact according to the result of this study, the necessary 
r>ser height Is mainly dependent on the active fuel lenpth and the 
core poufer density. 

The above results are much In agreement with the published 
results of Yasuso [13 and Kataoka t23 as 1t Is clear from 
(5. 1 ) and (5.2). 
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Table (5.1): Variation of decay ratio with coolant flow per bundle. 


flow per 
bundle 
( Kq/s ) 


decay ratio 


for q=42 KW/1 for q=50 KW/1 


results results In results results In 
obtained ref. Ill obtained ref.LI] 


.729166 


.819128 


.527683 


.617643 


.376329 


.466294 


.271967 


.361939 


.198905 


.288884 


. 146693 


.236676 


. 108727 


.198714 


.080748 


170738 


.059928 


.149920 


.044334 


.134328 


.032614 -03 


.122609 


.023799 


,113795 


.017180 


.107177 


.012234 


. 102232 


.008566 


.098565 















50 


Tsblft (5.7): Variation of mlnlnun cor® ■flow and necessary riser 


helglit with active fuel length. 


active 

fuel 

length 

( In m ) 

for c 

' 42 KW/1 

for q 

« 50 KW/1 


nilnlmm G 
per 

bundle 

' 

riser height (m) 

mlnlnxjm G 
per 

bundle 

riser height (m) 

results 

obtained 

results 
In ref[l 

results 

obtained 

results 
in ref Cl' 

7.3 

3.75 

2.4649 

- 

5.25 

2.7219 

- 

2. A 

5.32 

2.8248 

3.0 

6.82 

3.2549 

3.5 

2.5 

6.73 

3.2504 

3.7 

8.23 

3.9147 

4.5 

2.6 

8.04 

3.7494 

4.3 

9.54 

4.7184 

5.8 

2.7 

9.30 

4.3297 

5.2 

10.80 

5.6842 

7.0 

2.8 

10.56 

4.9995 

6.1 

12.06 

6.8309 

8.0 

2.9 

11,87 

5. 7675 

7.0 

13.37 

8.1784 

9.1 

3.0 

13.28 

6.6424 

7.9 

14.78 

9.7476 

10.2 

3.1 

14.85 

7.6334 

8.8 

16.35 

11,5603 

11.3 

3,7 

16.62 

8.7502 

- 

18.12 

13.6397 

- 

3.3 

18-66 

10,0026 

- 

20.16 ! 

16.0099 

- 


\ 
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active -fuel length (m) 

n,.<s.«.COOLRNT FLOW vs RCTIVE FUEL LENGTH 
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Table (5.3) Variation of riser height with core power density. 


power 

density 

(KW/1) 

riser height 1 

n ( meter ] 


for 

L*2.^. 

for 

L*2.7m, 

for 

L=3,0in. 

for 

L*3.3m, 

40.0 

2,7596 

4.2433 

6,5382 

9,8866 

40.5 

2,7655 

4,2506 

6,5475 

9.8989 

41.0 

2.7856 

4.2791 

6.5848 

9.9451 

41.5 

2.8053 

4,3050 

6.6152 

9.9772 

42.0 

2,8248 

4.3297 

6.6424 

10.0026 

42.5 

2.8443 

4.3544 

6.6696 

10.0284 

43.0 

2.8490 

4.3554 

6.6624 

10.0071 

43.5 

2.8692 

4.3841 

6.7005 

10.0550 

44.0 

2,8899 

4.4169 

6,7492 

10.1246 

44.5 

2,9115 

4.4550 

6,8122 

10.2231 

45.0 

2.9341 

4,5000 

6.8931 

10.3582 

45,5 

2.9580 

4,5532 

6.9958 

10,5377 

46.0 

2,9834 

4,6162 

7.1247 

10.7704 

46.5 

2,9949 

4,6645 

7,2441 

11.0086 

47.0 

3.0239 

4.7519 

7.4386 

11.3768 

47.5 

3.0551 

4.8546 

7.6746 

11.8310 

48.0 

3.0889 

4.9746 

7.9585 

12.3861 

48.5 

3.1254 

5.1146 

8.2981 

13.0591 

49.0 

3.1650 

5.2773 

8.7023 

13.8704 

49.5 

3.2080 

5,4660 

9.1814 

14.8444 

50.0 

3.2549 

5.6842 

9.7476 

16.0099 






RISER HEIGHT vs CORE POWER DENSITY 





CHAPTER 6 
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RECOWENOATION FOR FUTURE WORK 


1 . 


2 . 


3. 


4. 


5 . 


6 . 


In tTie stability analysis, the perturbation 1n the Inlet 
entltalpy of the coolant was neglected assuming that it will 
not make much effects on the results. But a more accurate 
results can be obt a ined by iivjluding inlet enthalapy 
perturbation in the analysis. 

The uniform axial heat flux distribution was considered in 
stability analysis. However the peaking factor that gives 
conservative results, was used but it will be desirable to use 
the actual distribution that is the sinusoidal distribution 
for this analysis. 

In calculating riser heij^t one dimensional model is 
considered, but a more accurate results can be obtained 
using two-dimensional model with actual heat flux 
distribution. 

In riser height calculation more accurate results can 
be obtained by considering pressure drop components of the 
bodies submerged in the core such as spacers, separators, 
plates, etc., which were not considered in present analysis. 
The analysis has been done assuming uniform distribution of 
coolant flow among bundles and than considering the worst 
bundle but it will .be better to consider the flow 
distribution corresponding to radial heat flux distribution 
so that stability criteria can be applied in each bundle. 


The effects of parameters like reactor pressure. 


number of 


bundles, core diameter and 
height which are not considered 


roaffiy others on riser 
in present thesis , can be 


studied in future. 
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APPENDIX A 

REACTOR SPECIFICATION 


Parameter 

Value 

Operating pressure (p) 

7.0 MPa. 

Active fuel length (L) 

2,3 to 3.3 m. 

Fuel rod diameter (d) 

12.3 mm. 

Core diameter (D) 

5.8 m. 

Volumetric power density (q) 

40 to 50 Kw/1. 

Core Inlet subcooling 

33 KJ/Kg. 

Fuel bundle type (M) 

8X8 square. 

Width of bundle (W) 

0.13813 m. 

Number of fuel bundle (N) 

872 

Local peaking factor (f^^) 

1.2 

Axial peaking factor (f^) 

1.5 

Radial peaking factor (f^^) 

1.4 

Void- reactivity coeffleclent (K^) 

- 0.08 % / X void. 

Flow rate per bundle (G) 

5 bo 20 Kg/s. 

Fuel material type 


Delayed neutron fraction (p) 

0.0066 

Decay constant of precursor (X) 

0.075 / sec. 

Promt-neutron life time (A) 

0.001 sec. 


Table A. 1 : Reactor Specification. 
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APPENDIX B 

PROGRAM FOR CALCULATING DECAY RATIO 

REAL L,MU,K,KV,LEMDA,MT,MG,NBH,NBH1,KO,MF,NU,NN(10),ND(10) 
REAL NTlN(10),NT1D(10),Mr2N(10),NT2D(10) 

D]>ENSION Y1N(10),Y2N(10),Y1D(10),Y2D(10),T1N(10),TTD(10) 
DIMENSION T2N(10),T2D(10),P1N(10),P1D(10),P2N(10),P2D(10) 
DIMENSION F 1N( 10) ,F1D( 10) ,F2N( 10) ,F2D( 10) .GEhCC 10) ,GEU6( 10) 
DIMENSION GEN(10),GEN1(10),(^N2(10),GEN3(10),GEN4(10) 
DIMENSION HN(10),HD(10)»Qa(10),GD(10),GiHN(10),GHD(10) 

{ CONSTANT DECAL ERATION } 

HF=1258.18 

HG=2769.67 

DHI=33 

VF=0.0014 

VG=0.028 

P=1000 

RELOSS=0. 1 

QV=42.E03 

O=12.3E-03 

L=2.7 

PF= 10970 

CF=171.68E-03 

MLN2.48E-03 

F=0-025 

PC=750 

G=9.81 

' BETA=0.0066 

GAMA=0.001 
LEMDA»0.075 
KV^O.OB 
W=0. 13813 
N=872 
M=64 
DC=5.8 
PIE=3.1414 
PR=0.35 
K=10.0 
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DO 5 GPB=1,20 
MT=GPB*N 
HFG=HG-HF 
VFG=VG-VF 

{ TO GET EXIT QUALITY OF STEAM } 

MG=TPG/ ( HGTRE LOSS ) 

XE=MG/MT 

{ TO GET NDN-BOILIN3 HEIGTH ) 

HI=HF-DHI 

NK^1=DHI/(DHI+XE*HFG) 
hBH= (L*ACOS( 1 ) )/PIE 

{ TO ^T COOIjANT flow VELOCITY IN NON-BOILING PART ) 
AP=PIE*D 

AX= (W*W-M*PIE*D*D/4 )/M 

AA=AP*L 

GO=MT/(M*N) 

V10=G0/AX 

QS=QV*D/4 

C0=(0.71+0.0001*P)-1 
KO=VF/(VGrCO*(XE+VF*( 1-XE )/VG)**2) 

{ TO GET HEAT TRANSFER COEFFICIENT ) 

CW=4 . 0* ( W*W-M*PIE*D*0/4 . 0 ) / ( 4 . 0*W+M*^IE*D ) 
RE=PCr*V10*D/MU 

M3= ( K*0 . 23*R£**0 . 8’^PR**0 . 3 ) / ( 1 0O0’*'D ) 

«>QS*AP*VFG/ ( FFGTAX ) 

TEXIT=- (ALOG(V 10/(Lr(L-NBH)+V 1 0 ) ) ) 

MF=PIE'*D*0*L’W/4 
MJ=PC*AX*DHI/ (QSf^AP ) 

B0=BETA/ ( LEMDA*GAMA ) 

{ TO GET Ms) } 

Y1N(1)=-N8FrH0*AA 
Y 1 N (2 ) =-0 . 2*NBH*MF*CF 

Ijirv / t \ i A A 
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GEN2(3)=1 
CALL IMT(GEN2,3) 

GEN3(1)=1 
(^N3(2)=-NU 
CALL INT((2N3,2) 

GEN4(1)«E 

GEN4(2)=CC 

CALL INT(GEN4,2) 

CALL MULTYCGENI . GENS, GENS) 
CALL MDLTY(GEM4,GEN2,GEN6) 
CALL ADD(GEN5,G£N6,(SN1) 
CALL PROO(Y,Y1N,(2N4) 

CALL MULTY(GEN4,GEN2,GEN) 
CALL MJLTYCGENl ,Y1D,GEN4) 
CALL AD0(GEN,GEN4,P1N) 
CALL MLJLTY(Y1D,GEN2,P1D) 

{ TO GET NTCs) } 

NN(1)=BO 
fW(2)=B0 
CALL INT(NN,2) 

ND(1)=0 

ND(2)=BO+1 

ND(3)=1 

CALL INT(ND,3) 

{ TO GET H(s) } 

CALL ADO(F2N,P2N,GENl) 
GEN2»F2D 

CALL MJLTY(F1N,P1D,GENS) 
CALL MLJLTY(F1D,P1N.GEN6) 
CALL ADO(GEN5,GEN6,GEN3) 
CALL MJLTY(F1D,P1D,GEN4) 
CALL MJLTYCGENl ,GEN4,HN) 
CALL MJLTY(GEN2,GEN3,HD) 


{ TO GET G(S) ) 



CALL WJLTY(hKD,T2D,NT2D) 

CALL PR0D(KV,m'2N,GEN) 

CALL SUB(NT2D,GEN,GEN1) 

CALL MJLTY(NTlN,hn-2D,GEN2) 

CALL PR0D(KV,G£N2,GN) 

CALL MULTY(GEN1,NT1D,GD) 

{ TO GET G‘H(s) } 

CALL ML)LTY(GN,LW,GHN) 

CALL MJLTY((a),HD,aC) 

{ TO GET DECAY RATIO } 

A=GD(7)/GN(6) 

B= {GD(6 )-A*GN(5) )/GN(6 ) 

POLE1=-B/A 

A=HD(7)/HN(5) 

B= ( - HD ( 6 ) + A*HN ( 4 ) ) /HN ( 5 ) 

0(HD(5)-A*HN(4)-Bi*HN(3))A«(5) 

D0=B*B-4*A*C 

ALPHA=-B/(2*A) 

BETA1=SORT(-DO)/(2*A) 

TP=(2^IE)/BETA1 

DR=EXP(-ALPHA'nP) 

PRINr*=,GPB,DR 

5 CONTINUE 
STOP 
END 

{ SUBROIJnNE TO INITIALIZE A POLYNOMIAL } 
SUBROLTFINE INT(POLY,N) 

DIMENSION POLY (10) 

DO 60 I=N+1 ,10 
POLY(I)=0.0 

60 CONTINUE 
RETURN 
END 

{ SUBROUTINE TO MULTIPLY TWO POLYNOMIALS } 


n-Tk*- fcji fl TV Z' on* V 1 DTi V9 DW tt T ^ 
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DO 80 1=1 ,10 
PMULTT(I)=0.0 
DO 70 J=1,I 

PMJL TT ( I ) = PWJLTT ( I ) + POl Y 1 ( J )*Pa Y 2 ( I - J + 1 ) 

70 CONTIMUE 

PMLJLT(I)=PHJLTT(I) 

80 CtWINLlE* 

RETURN 

END 

{ SUBROUTIIC FOR PRODUCT OF POLYNOMIAL WITH A SCALAR } 
SUBROUHNE PROD (X, POLY , PRO) 

DIl^NSION Pa_Y(10),PRO(10) 

DO 90 1=1,10 
PRO(I)=POLY(I)’^X 
90 CONTINUE 
RETURN 
END 

{ SIBROUTINE TO ADO TWO POLYNOMIALS } 

SUBROUnNE ADO(POLYl,POLY2,AD) 

DIMENSION POLY1(10),POLY2(!0),AD(10) 

DO 95 1=1,10 

AD( I) =POL Y 1 ( I ) +POL Y2 ( I ) 

95 CONTINUE 
RETURN 
EhD 

{ SUBROUHNE TO SUBSTRUCT TVO POLYNOMIALS } 

SUBROUTINE SUBCPOLYl ,POLY2,SU) 

DIMENSION POLY1(10),POLY2(10),SU(10) 

DO 99 1=1,10 
SU(I)=POLYl(I)-POLY2(I) 

99 CONTINUE 
RETURN 
END 
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APPENDIX C 

SUBROUTINE FOR CALCULATING RISER HEIGHT 

{ SUBROUTINE FOR CALCU^TING RISER HEIGHT FOR GIVEN FLOW PER ) 

{ BUNDLE, ACTIVE FUEL LENGTH, AND CORE POWER DENSITY } 

SLBROUTINE RISER(GPB,L,QV,RH) 

REAL L,M,N,MT,MG.NBH1,>©H,KR 

W^=1258. 18 

HG=2769.67 

DHI=33-0 

VF=O.OOU 

VG=0.028 

R£LOSS=0.1 

D=12.3E-03 

F=0.025 

PC=750.0 

G=9.81 

W=0. 13813 

N=872 

M=64 

PIE=3.1415 

S=1.9 

FR=0.025 

F0=0.025 

D05.8 

KR=0.5 

TPG=QV* ( PIE*DC*DC*L/4 . 0 ) 

MT=(yB*N 

HFG=HG-HF 

VFG=VG-VF 

eT=TPG/MT 

{ TO GET EXIT STEAM GUALITY } 

MG=TPG/(HGr ( 1 -RELOSS) ) 

XE=MG/MT 

{ TO GET NON-BOILING HEIGHT } 

HI=HF-DHI 



NBH= (L*ACOS( 1 .0-2. 0*NBH1 ) )A>IE 

{ TO GET FLOW VELOCITY IN NON-BOILING PART OF CHAWEL } 
AX= (W*V-M*PIE*D*D/4. 0 )/M 
G0=MT/(M*N) 

V10=G0/(AX*PC) 

{ TO GET HYDRAULIC DIAMETER OF COOLANT FLOW } 
DO=DC-‘SORT(1.0+KR) 

DH=4.0*(W^i;-^rPIE’^D/4-0)/(4.0*WrPIE*D) 

{ TO GET FRICTIONAL PRESSURE DROP IN CHAMnEL } 
SIE=S*(VF/VG) 

AE=XE/(XE+(1-XE)’^IE) 

R= ( 1 + ( 1 -AE )** (- 1 )+ ( 1 -AE )*=* (-2 ) )/3 . 0 
DPF=F*V 1 O^V 1 0*PC* (NBH+R* (L-NBH) )/(2 . 0*DH) 

{ TO GET ACCELERATION PRESSURE DROP IN CHANNEL } 

X= ( ( ( 1 -XE )**2 )/ ( 1 - AE ) +XE*XE*VG/( AE*VF )- 1 . 0 ) 
DPA=PC*V10*V10*X 

{ TO (^T PRESSURE DROP AT ABRUPT AREA CHANGES } 

A2=N* (Vnrf-«*PIE*D*D/4 . 0 ) 

A1=PIE*DC*DC/4.0 

DPC=0. 7*PC^V10*VlCr( 1-(A2/A1 )**2) 

DPE=Mr»MT*VF*(X+1 )*( (A1*A2 )**'(-! )-A1**(-2)) 

{ TO GET AVERAGE DENSITY IN THE BOILING PART OF CHANNEL } 
Cl = (QT-2 . 0*DHI )/(2 . 0*H^G) 

C2=-QT/(2.0*HFG) 

C3=SIE+(1-SIE)*C1 

C4=C1-SIE)’*C2 

pe=aeag+(i-ae)/vf 

PCE=PC-PE 

PG=VGr*(-l) 

IF (C3=**2.GT.C4**2) THEN 
X=PIE*NBH/(2.0*L) 

X= ( (C3-C4 )’^IN(X )/COS(X ) )/SGRT (C3**2-C4**2 ) 

X= 1 . 0-2 . O^ATANC X ) /PIE 

X=(C1’^:4-C2*C3)*L*X/(C4*(L-NBH)*SQRT(C3**2-C4**2)) 



ENDIF 

IF (CA**2.GT-C3**2) THEN 
X=PIF*NBH/(2.0*L) 

X= (C^-C3)*SIN(X )/COS(X) 

X= ( X+SQRT (C4**2-C3**2 ) )/ ( X-SQRT(C4**2-C3**2 ) ) 

X= (C 1*C4-C2*C3 )*L’*'ALOG(X )/(C4*(L-NBH)*PIE*9QRT (C4’**2-C3**2 ) ) 
PCB= ( PC-PG )* (C2/C4-X ) 

ENOIF 

{ TO GET VELOCITY OF COOLANT IN RISER ) 

VE = 4 . 0*MT / ( PE’^PI E'^DC^DC ) 

{ TO GET VELOCITY OF COOLANT IN DOWNCOMER } 

VO=4 . 0’TflT/(PC*PIE*(D0*D0-DC^) ) 

{ TO GET RISER HEIGHT } 

A1 1 =FR*‘VE*VE*PE/( 2 . 0*IX )+FD*VD*V[m:/(2 . 0* (DO-DC) ) 
Cn=Fiy*'VD*VO*F<rL/(2 .0*(D0-DC) ) 

D1=PCB’*'(L-NBH)*G 

B1=C1 1+Dl+DPF+DPA+DPC+DPE 

RH=B V(PCE*G-An ) 


RETURN 
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